Scaling laws for precision in quantum interferometry and bifurcation landscape of optimal state 
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Phase precision in optimal 2-channel quantum interferometry is studied in the limit of large photon number 
N ^ 1, for losses occurring in either one or both channels. For losses in one channel an optimal state undergoes 
an intriguing sequence of local bifurcations as the losses or the number of photons increase. We further show that 
fixing the loss paramater determines a scale for quantum metrology - a crossover value of the photon number 
Nc beyond which the supra-classical precision is progressively lost. For large losses the optimal state also has a 
different structure from those considered previously. 

PACS numbers: 42.50.-p,42.50.St,06.20.Dk 



It has been recognized that using quantum states of light 
may increase the resolution of interferometric measurements 
lUKH. Particular states of N photons achieve the Heisenberg 
limit of phase resolution for standard error on the phase esti- 
mate Ai^ = 1 /N, an improvement over the classical (or shot- 
noise) limit Aip = 1/VN that is obtainable when N photons 
enter the interferometer one at a time. These bounds are de- 
rived by an application of the Cramer-Rao inequality |2 | for 
the standard error of an unbiased estimator, Aip > [vF)^^/'^, 
where F is the quantum Fisher information (QFI) |4| and v 
is the number of repeated independent trials. Assuming any 
instrument is composed of three components: quantum input 
state, dynamics and measurement; the functional F depends 
only on the first two — it assumes an optimal measurement 
choice. For pure states in a single mode F/A = AP"h = 
{n^) — (fi)'^ (where h is the number operator) and a famil- 
iar uncertainty relation is recovered: AnAtf > 1/2. Thus, 
for a lossless two-mode interferometer QFI and precision are 
greatest for the maximum variance state, or 'NOON state' ; it 
saturates the Heisenberg limit. Unfortunately, it is also highly 
susceptible to noise, especially dissipation |5]. 

To mitigate this problem various two-component states 
were proposed f6'-'8l, where the loss of a number of photons in 
the first mode does not destroy the superposition. The preci- 
sion performance under dissipation of various Gaussian states, 
e.g. squeezed, coherent and thermal states, has also been con- 
sidered recently |i9J. In all cases, the precision was found to 
be supra-classical for certain range of losses and N. 

In the lossy case the pure input state of two oscillator modes 
maximizing QFI 

N 

10) = ^0„|n)i|7V-n)2, (1) 

n=0 

must balance supra-classical precision against robustness to 
photon loss. In this notation the NOON state has two non-zero 
components, 4>q = cf)iq = l/y/2. For a lossy interferometer 
light propagates in each arm as a damped harmonic oscillator, 
with frequencies w*^^-',a;'^^ and dissipation 7*^^'', 7^^^. Equiv- 
alently, losses can be introduced by beam-splitters in each 
mode with reflectivity i?^^^^) = 1 - exp{-7(^'^'t}. Those 
lost photons siphoned out of the modes are then traced over. 



In the simpler case of losses in only one of the two modes, 

as might occur when that mode is 
directed through a partially transparent test sample, the state 
10) decays into a mixture p = |V'fc)(V'fc| with 

l^fc) = ^Ev^e™'^'/'"!"-^'^-^)' (2) 

corresponding to the loss of k photons. Here Wk is the nor- 
malization factor; the phase differense is ip — (cj^^^ — a;^^^) t 
and the loss enters via coefficients A„:fc — Ry^^'^. 

Fisher information of the mixed state resulting from losses is 
a weighted sum over pure components ^ = X^fc WkFk, where 
Fk = 4A^ni for pure states jV'fc) (6\. 

Refs. |6, 8 1 use numerical optimization to construct states 
|0opt) maximizing this information for a range of photon 
numbers N and loss R. For larger photon numbers/loss, 
a complicated multicomponent structure arises, but the evo- 
lution of the optimal state with increasing loss was previ- 
ously not well-understood. More importantly, the question 
of the asymptotic scaling with large N of the Fisher infor- 
mation has been left unanswered. The significance of this 
issue becomes apparent when one examines the suitability 
of quantum-enhanced sensors for tasks such as gravitational 
wave observation where ^ 1 is necessary to reach the de- 
sired sensitivity. Later we will discuss the case of free-space 
target acquisition and ranging (e.g. quantum LIDAR) where 
combined loss R due to atmospheric attenuation and limited 
target reflectance is typically over 99%. 

Here we study the analytically tractable limit ^ 1 by 
treating n/N = x £ [0; 1] as a continuous parameter. Exam- 
ining the limits of small and large loss has revealed a scaling 
relationship for the optimal Fisher information 

F,,,{N,R)=N^F(^^^y (3) 

that cleanly interpolates between these limits. The non-trivial 
dependence on A^ and R is captured by a single quantity: r = 
NR/{1 — R). The structure of the optimal state also depends 
on r alone, save for small differences due to discrete nature of 
parameter x = n/N. We were able to demonstrate that for any 
finite r the optimal state can have only a finite number of com- 
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ponents. This number increases with r as the optimal state un- 
dergoes a sequence of bifurcations: unbalanced NOON state 
ceases to be optimal for r < r[ Ki 0.912957 QO], superseded 
by a state ^r^|iV)i|0)2 + y^\xiN)i\il - Xi)N)2 as 
has been noted in f6l (a similar state has been proposed in 
Ref |7|). For larger values of the paramater (r > r2), an op- 
timal state acquires a third component |0)i|A^)2, which shifts 
away from the origin to |a;2A^)i|(l — X2)N)2 for r > r'2. The 
universal set of bifurcation points r'^ < r2 < r'2 < r^^ < r'^ < 
• • • as well the weights and locations xi of components in 
m + 1-component state \(f)) = ^~pi\xiN) are determined 
by solving a system of 2m — 1 or 2m equations. The re- 
sults are shown in Fig. [T] An important caveat is that since 
Ivh = xiN are not integers in general, single non-integer 
components split into two adjacent integer components for fi- 
nite N. 




FIG. 1: Probability weights p^, represented as stacked histograms 
(top) and locations xt (bottom), as a function of r. Different compo- 
nents are indicated using color. Black solid line on the bottom figure 
is the reseated Fisher information J-(r) and the blacked dashed is its 
asymptote 4/r; convergence takes place for much larger values of r. 
Thresholds r2 , ra , r4 correspond to appearances of new components 
at the origin. Components separate from the origin at critical values 
f''i,'i'2,r'o,,r^. The data for 20 ^ r ^ 25 is magnified (top figure, 
lower right corner) showing components with very small weight. 

The numerical results of Ref. fSl correspond to the leftmost 
region r ~ 1 of Fig. [T] In this work we are primarily inter- 
ested in the regime r \ (not shown). Our motivation is that 
the loss parameter R is determined by the properties of the 
medium and is, therefore, fixed. With R — Q increasing N 
provides continuous improvement in phase precision per pho- 
ton - we are interested in whether this remains the case in the 
presence of losses. (It should be noted that no efficient pro- 
cedure for generating optimal iV-photon states for very large 
N presently exists.) Many applications, such as gravitational 
wave detectors, pursue orders of magnitude improvement over 



current precision limits, hence the importance of finding out if 
any fundamental limits to the strategy of increasing the photon 
number N exist. 

With increasing N, optimal states seemingly increase in 
complexity as the value of r and the number of components 
increases. However, as the density of these components in- 
creases, for large N the optimal state |</)) may be approx- 
imated by a continuous function. Examination of Eq. ([3]) 
shows that, r being proprtional to N, a quadratic (Heisenberg) 
scaling of the Fisher information is obtained only for small r, 
where a{r) — — d(ln J^)/d(lnr) is small, see Fig.[l] (Here 
we are fitting to a local exponential model: F ^ r^".) In 
the opposite limit r ^ 1 one has F ^ 4/r, i.e. a = 1. so 
that for iV ^ 1 the Fisher information increases only linearly, 
F oc N. Precision, quantified by F, will always improve 
at least linearly with the photon resource N, (by sending the 
photons through the instrument one at a time) but the more in- 
sightful question is how does the amount of 'intrinsic' Fisher 
information, i.e. per photon, scale with N if photons are com- 
bined in some optimal quantum superposition - indeed, how 
does this scaling change, given a loss parameter i?? 

General Upper Bound: This general linear upper bound can 
be demonstrated without making approximation of large N or 
R ^ 1. Since the variance is unaffected by a constant shift, 
one can rewrite the F as 



R 



(4) 

Observing that the second term is negative and performing the 
sum over k in the first term, we obtain the inequality 



R 



R 



(5) 



This upper bound (J^ippei) is always valid for any R > 
demonstrating that quadratic precision (at the Heisenberg 
limit F oc N'^) is only possible for R oc When R is 

fixed, it implies T oc N, scaling proportional to the shot-noise 
limit. This bound also appeared recently in the complemen- 
tary context of global phase estimation 1 11 1. 

Limit J-"upper is reachable asymptotically as can be shown 
by constructing a wavefunction that minimizes the correction 



A J- = J-,,, 



T. In the limit i? ^ 1 we approximate the 



true optimal state {(j>n} by a continuous function cl){x) (with 
X = l — x = 1 — n /N), smooth on scales ^ 1 / \/N and obtain 
approximately 



4Ar2 .00 
Jo 



r 



dx. 



(6) 



where the upper limit has been set to infinity since the width of 
(l){x) is much smaller than 1. The term proportional to (f>^{x) 
is the first term of Eq. Q subtracted from J-"upper, and the term 
proportional to 4>''^{x) is the second term in Eq. ^ taken with 
the opposite sign. 
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Minimization of (|6]l subject to the boundary condition 
4>{0) — 112 | and the normaHzation constraint produces 



(r/4) 



1/6 



Ai'(Mi) 



Ai 



1/3 



Ail 



(7) 



where Ai(2;) is the Airy function, /ii w —2.338107. . . is its 
first (largest) zero, and the prefactor ensures normaUzation. 
Together with the next order correction, the Fisher informa- 



tion of the optimal state Jopt = J"i 
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For i? - 1 (r - N) the width of ^ is 0{N'^/^); so is the 
leading correction in Eq. (|8]l. The upper bound becomes sat- 
urated when the number of photons exceeds a value of Nc 
estimated by equating the principal term and the leading order 
correction in Eq. ([s]). This yields Nc — rc(l — R)/R with 
r-c-4|^i|3~50. 

Arbitrary Loss in Both Arms: Whenever both i?*^^^ and R^^^^ 
are non-zero, the density matrix is a mixture of pure states 
\1pk1k2) resulting from the loss of fci and fc2 photons in modes 
(1) and (2) respectively. In Eq. (|2]) the factor y^Aj^ becomes 

\J ^n-ki^n-k2 stateslV'fe) become liAfcifca)- 

The number of photons lost in each mode is not observed 
directly, although their sum k = ki + k2 can be inferred 
by subtracting the detected photon number from the input N. 
Consequently the linear decomposition of the Fisher informa- 
tion serves only as an upper bound J' ^ J2ki k2 "^kik^-^kik^ 
lUl and the determination of quantum Fisher information re- 
quires the diagonalization of the density matrix l?): 



F ^ AY,\,{v.,\nl\v,) - 



E 



8A2 

Ai + \j 



\{v,\n^\v,)\\ (9) 



where Ai and \vi) are eigenvalues and eigenvectors of the den- 
sity matrix respectively. Diagonalizations within subspaces 
corresponding to a fixed total number of lost photons k = 

ki + k2 may be carried out independently. In the limit N ^ 1 

(1 2) 

the coefficients AJ^.j, may be approximated by Gaussians so 
that the corresponding density matrix is also Gaussian in the 
continuous limit as long as i?*^^' ^ i?^^) and the wavefunction 
is smooth on scales ^ y/N. This density matrix may be 
expanded in terms wavefunctions of harmonic oscillator with 
the aid of Mehler formula |13 | and the sum (|9]) is evaluated 
noting that non-zero matrix elements correspond to j = i ± 1. 
The surprising outcome is that the exact Fisher information 
equals the linear upper bound ( J^ippei for arbitrary loss in both 
arms) in the asymptotic limit. This is also true in the sym- 
metric loss case (i?*^^^ = i?^^^) as the optimal state itself turns 
out to be a Gaussian. This case has some import; firstly, it 
is relevant for balanced instruments where phases may be in- 
troduced in either arm, e.g. gyroscopes, and secondly; the 
analysis has an extended applicability beyond losses in modes 



(1) and (2) to those occurring in any superposition of these 
modes. Accordingly, the discussion is applicable to losses in 
detection after the mode mixing. 

Expressed in terms of parameters r^^'^") = 
NR^^'^y (l - R^'^'^^), the upper bound Q changes to 

The optimal wavefunction is computed by minimizing the 
correction to the Fisher information [x = n/N, x^, — 




dec, 



(10) 



which produces a Gaussian centered at x = of width 

y2(r(i)r(2))i/8/(y^ + Vr(2y). This width scales as N^/'^ 

(Cf Ar2/3 for 

single mode losses). 
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FIG. 2: Optimal 20-photon states for 95% loss (r — 380) in one 
(left) and two (right) modes. Red bars represent amplitudes </!>„ ob- 
tained by numerical optimization. Black lines represent analytical 
approximation with an Airy function and a Gaussian. These optimal 
states offer a precision improvement (square root of Fisher informa- 
tion) over coherent light of just 6% (single-mode losses) and 0.4% 
(symmetric losses), owing to high loss amount. (For single-mode 
losses, using coherent light, the precision used in calculation is for 
the optimal reflectivity of the beamsplitter in MZ interferometer.) 

For moderate losses, this optimal form is attained when the 
number of photons is large. In the limit of large losses, this 
asymptotic form is reached with a small number of photons 
(see Fig. |2]i, within reach of current laboratory capabilities. 
The Fisher information together with the leading correction is 




opt 



The correction scales as A^^/^, in contrast to the A^^/^ scal- 
ing for single mode losses. CoiTespondingly, the crossover to 
the limiting behavior is expected for smaller A^. The conver- 
gence to asymptotic precision for the case of single mode and 
symmetric losses is illustrated in Fig. [3] 

Summary and Outlook: An important question in interfer- 
ometry is that of scaling of precision with photon number N, 
and whether this scaling changes as N increases. We have 
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FIG. 3: Left: Fisher information for symmetric (ii*^^' — i?'^' = R) 
and single-mode (i?'^' = R, i?'^' — 0) losses for optimal A'^-photon 
states as a fraction of the linear upper bound [N{1 — R)/R and 
4A'^(1 — R)/R respectively]. The curves must tend to 1 for large 
A'^, but the convergence is faster for losses in both arms. For sym- 
metric losses we use exact Fisher information, not the approximate 
upper bound. Right: Collapse of data when replotted as a function of 
r = NR/{1 - R) [see Eq. ^] 



shown analytically that the best supra-classical precision, as 
quantified by the QFI, is quadratic in N initially and under- 
goes the crossover into the regime where it scales linearly with 
N. This has far-reaching repercussions — the crossover loca- 
tion Nc is a function only of dissipation, due to these losses a 
resource bound (or scale) is imposed on any instrument claim- 
ing to offer supra-classical precision. Our results, indicat- 
ing that the improvement of precision is dissipation-limited, 
should temper expectations of various proposals advertising a 
quantum leap in sensitivity of gravita- tional wave detectors 
with non-classical light 1 15l [T6l . 

We obtained the scaling relationship for the Fisher infor- 
mation and studied the evolution of the optimal state: it un- 
dergoes a sequence of bifurcations for intermediate values of 
N. In the limit of very large N, the optimal state becomes a 
continuous function with width that scales as N^/^ or N^^^ 
(for one/two mode loss), a signature property that makes it 
distinct from states considered previously. 

It should be noted that the precise value of the exponent e 
for the width scaling (or the exact form of the wavefunction) 
does not affect the asymptotic value of the Fisher information. 
It is vital, however, that the width be much greater than y/N 
and much smaller than N. Any exponents lying strictly within 
1/2 < e < 1 will give the asymptotic value of the Fisher 
information, though convergence rates will be optimized with 
2/3 or 3/4. The spin coherent state discussed in Ref. llT4l 
with respect to the large loss limit has a width scaling cx N^^'^ 
and is suboptimal - it does not offer supra-classical precision 
for any losses in view of Caves' theorem |15|. Similarly, the 
Holland-Burnett state ifTTl with width oc iV is also suboptimal: 
for symmetric losses its asymptotic Fisher information is half 
the upper bound. 

As a final note, we enlarge upon two proposed applica- 
tions of quantum light; to free-space target acquisition and 
ranging, and to gravity wave observation. In clear weather. 



infra-red light is attenuated by 0.5 — l.OdB/km. Locating a 
target at 10km distance requires 20km of roundtrip propaga- 
tion, i.e. i? — 90 — 99% loss, combined with a typical 10% 
target reflectance becomes 99 — 99.9%. Optimal Fisher infor- 
mation per received photon is 4/i?, thus naively one would 
expect a two-fold improvement in phase precision 5(j) over 
coherent light (having Fisher Information unity per received 
photon) even for high losses. This comparison is for an in- 
terferometer with 50/50 beamsplitters; precision can be triv- 
ially increased with coherent light inputs by optimizing the 
beamsplitter reflectances. Compared with this strategy, non- 
classical light can improve precision by at most a factor of 
(1 + Vl^ni)/VR, i.e. by 3 - 10% for losses above. This 
fractional advantage in the very high loss limit does not offset 
the high practical cost of generating those optimal states we 
have discovered. This result should moderate expected out- 
comes of such proposals. 

To contrast, consider two-mode losses R w 1%, the ex- 
pected domain of advanced interferometric gravitational wave 
detectors with high-reflectivity miiTors and state-of-the-art 
photodetectors. The improvement to Scf) over classical light 
for the same ^ 1 approaches a factor of This 10- 

fold improvement falls far short of more optimistic estimates 
assuming idealized conditions 1 16 1 but still represents a clear, 
non-trivial advantage for the optimal input states we have dis- 
covered. The ability to reduce intrinsic quantum noise by an 
order of magnitude without an associated increase in radiation 
pressure noise (the photon flux has not increased) is certainly 
of interest for the development of gravity wave detectors. 

We remark that due to the isomorphisms between two har- 
monic oscillators and spin/qubit systems these results apply 
quite broadly in metrology protocols, from photonic systems 
and atomic condensates, to spin ensembles coupled to heat 
baths, and other processes undergoing both unitary evolution 
and dissipation. 

G.A.D. contributed to this work while under contract with 
Mission Critical Technologies, Inc. 
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